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Abstract 



The Kazakov-Migdal (KM) Model is a U(N) Lattice Gauge Theory with a 
' Scalar Field in the adjoint representation but with no kinetic term for the 

cn 

' Gauge Field. This model is formally soluble in the limit N ^ oo though ex- 

o 

l/^ ' plicit solutions are available for a very limited number of scalar potentials. A 



"Double Penner" Model in which the potential has two logarithmic singulari- 



Q^' ties provides an example of a explicitly soluble model. We begin by reviewing 

' the formal solution to this Double Penner KM Model. We pay special at- 

^ - tention to the relationship of this model to an ordinary (one) matrix model 

' whose potential has two logarithmic singularities (the Double Penner Model) . 

We present a detailed analysis of the large N behavior of this Double Penner 
Model. We describe the various one cut and two cut solutions and we discuss 
cases in which "eigenvalue condensation" occurs at the singular points of the 
potential. We then describe the consequences of our study for the KM Model 
described above. We present the phase diagram of the model and describe its 
critical regions. 
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I. INTRODUCTION 



Several years ago Kazakov and Migdal proposed a model which they hoped would 
provide a description of Quantum Chromodynamics (QCD) in the limit in which the number 

of colors is large. Their model is a Lattice Field Theory defined on a D dimensional 
hypercubic lattice with sites labeled as x and links labeled by pairs of nearest neighbor sites 
< x,y >. The fields in the model are a scalar field which, for every site x, is an x 
Hermetian matrix and a Gauge field U{x, y) which, for every link < x, y > of the lattice, is 
an X Unitary matrix. The Action for a given configuration and f/(x, y) in their 
model is given by 

S^m = NY.1:iV[<P{x)]-N Tr(</)(x)[/(x,y)0(y)?7t(x,y)) (1) 

X <x,y> 

where V{(j)) is a potential function for the scalars which is, at this stage, arbitrary. The 
factors of A^ are included to assure a smooth large A^ limit for the model. The model is 
defined by the Partition Function Zkm given by the following Functional Integral: 

= / n ^^(^) n y)] [-'S'km] (2) 

X <x,y> 

Here d(l){x) is the Hermetian integration measure over the matrix and [dU{x, y)] is the 
invariant Haar measure for integration of the matrix U{x,y) over the unitary group U{N). 
This model is invariant under the gauge transformations 

(f){x) uj{x)(f){x)uj\x); U{xy) ^ uj{x)U{xy)uj\y) (3) 

where u{x) is an arbitrary U{N) valued function of x. The second term in Eq. (|lD is the 
usual gauge invariant kinetic term for a scalar field in the adjoint representation of a gauge 
group. In fact this model (which is called the Kazakov-Midgdal Model or the KM Model) is 
simply a model of an adjoint scalar coupled to a Gauge Field except that the usual kinetic 
term for the Gauge Field is omited. It is the absence of this term that makes the model 
soluble in the large A^ limit. We begin by reviewing how this works. 
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The first step, which can be done for finite A^, is to explicitly integrate over all the 
unitary matrices U{x, y) in Eq. (^. This can be done since there is no coupling between the 
various U{x,y). The result is a Functional Integral over the fields (f){x) which. As a result 
of the Gauge Invariance the integral depends only on the eigenvalues (f)i{x) {i = 1..N) 
of the matrices In fact the integral can be done explicitly using the "Itzykson-Zuber" 

formula 0. (For details see Ref. fl]].) The result is: 



Zkm oc / Yld(f)i{x)A'^[(f){x)] exp 



-iV^Try[</)(x)] 



n ^^^'^ ^ (A\ 



where d(j)i{x) is the ordinary integration measure over the real numbers and A[</)] = 

detij{(t>iy^^ = Y[i<j{(t^i ~ 4>j) is the Vandermonde determinant for 0. 

In this form of the Functional Integral it is well known [Q] how to go to the limit of 
large A^. When N oo the partition function is dominated by the stationary points 
of the action (defined as minus the logarithm of the integrand). In such circumstances the 
stationary points ($j(x)) are called the Master Fields for the Theory. The value of the 
integral when A^ ^ oo is equal to the integrand evaluated when 0j(x) is equal to the Master 
Field $i(x). 

When A^ is infinite it is customary and convenient to describe the Master Field will will 
usually be independent of a; by a density of eigenvalues. The idea is to order the eigenvalues 
so that <l>j is monotonically increasing and then to define the density of eigenvalues 

so that 

/oo 
pi\)d\ = 1 (6) 
-oo 

Np{X)dX is equal to the number of eigenvalues in a range dX about A. 

The simplest KM Model in which the potential V{(j)) = m^cj)^ /2 is quadratic was first 
solved by Gross 0]. There is a wealth of literature both on the solutions to the saddle point 
equations and on the relationship of the KM Model with QCD. A selection of references are 
Refs. i - [H . 



It was originally thought that this model would have a second order phase transition 
and some evidence was given in the case of a quadratic potential that this occurs when 

= 2D. It was then argued that the critical behavior of this model should be represented 
by QCD, the only known nontrivial four dimensional field theory with non-Abelian gauge 
symmetry. Unfortunately the solution of Gross ||^ showed that the Gaussian model had no 
critical behavior. This problem, combined with a further problem of an additional local 
symmetry which implies the vanishing of Wilson Loops and a better understanding of 
why the Gaussian Model fails to induce QCD has led to a consensus that the KM Model 
does not induce QCD. 

Nonetheless the model (Q) is interesting in its own right both as a Gauge Theory which 
is soluble in large and, as we shall review below, as an interesting example of a Matrix 
Model. Fortunately there has recently |21]| been some progress in finding an explicit 



solution to a non-Gaussian Kazakov-Migdal Model with a logarithmic potential of a very 
specific form. This solution was found by relating the density of eigenvalues for the KM 
Model to that of the an ordinary (one) Matrix Model whose potential has two logarithmic 
singularities. Although much is known about Matrix Models for polynomial potentials and 
for potentials with a single logarithmic singularity (the Penner Model) little is know about 
the model with two logarithmic singularities (which might be called a "Double Penner" 
Model). 

The purpose of this paper is twofold. First of all to study ordinary Matrix Models with a 
"Double Penner" Potential (which has two logarithmic singularities) and then to relate the 
solution of this problem to the solution of the Kazakov Migdal Model with a specific class 
of logarithmic potentials. 

The results of this paper are applicable in a much wider context than just the KM 
Model. If, for example, we choose the dimension D = 1/2 we recover the solution of the 
ordinary "Two Matrix Model" The techniques discussed in this paper and in Ref. |20 



are applicable to ordinary Matrix Models whose wide range of applicability can be seen in 
Refs. HH and m . 



The plan of the paper is as follows. In Section || we review the method of of finding 
the eigenvalue distributions for both Kazakov-Migdal Models and ordinary Matrix Models. 
We then show in general how this allows for explicit solutions in the case of a 'Double 
Penner" potential and how the density of eigenvalues for the KM Model is related to that 
of an ordinary matrix model. In Section |T| we describe in detail the variety of solutions 
to the ordinary Matrix Model with a Double Penner potential. We discuss the various one 
cut and two cut solutions and the cases in which there is "condensation of eigenvalues" at 



the singularities. In Section |^ we apply the results of Section |T| to the Double Penner 
KM Model. We present a phase diagram for the model and review its critical behavior. In 
Section we summarize our results and conclusions. 

II. METHOD OF SOLUTION 

A. Ordinary Matrix Model 

Consider first an ordinary Hermitian Matrix Model whose only variable is a single N x N 
Hermetian Matrix 0. The model is defined by by the partition function 



(f) exp[-NtiV{(t>)] (7) 

where d(j) is the Hermetian integration measure. If (as is always assumed) TtV{(I)) is invari- 
ant under U{N) transformations of then the integrand in Eq. (|^) depends only on the 
eigenvalues (pi of and the partition function can be written as 

Z (X Jlld(f)i A\(f)) exp[-NtTV{(f))] (8) 

i 

As discussed in the Introduction, in the limit N ^ oo this integral is dominated by a matrix 
$ whose eigenvalues /i are distributed according to some distribution p(yu). 

There are many methods of studying and solving such Matrix Models. The method 



we discuss here is most suitable for generalization to the KM Model []TT|. The density of 
eigenvalues can be found by defining the quantity 



where A is an arbitrary complex number. Once Ex has been computed the density of 
eigenvalues can easily be determined since at large N 

1 1 1 

N X — <P j-oo A — ?7 

So that 

27Ttp{r]) = Ex-ie-Ex+ie (11) 

which is nonzero only along the branch cuts of Ex. Thus by finding the branch cuts of Ex and 
computing the discontinuities across these cuts we can compute the density of eigenvalues 
p(A). 



Following Ref. |20| we begin with the equation 



This is now written in terms of an integral over the matrices rather than in tems of their 
eigenvalues. Eq. ( ]T^ ) leads to the following equation for Ex 

- (>''X^l"(^)> = (13) 

In the simplest case in which ^(</)) = m^(/)^/2 this equation leads immediately to a quadratic 
equation for Ex 

Ex = ]^ (m^A - Vm^A^ - Am?) (14) 
and thus to the well known semicircle distribution of the eigenvalues of 



p(A) ^ - ^ (15) 



which has support on the interval (— 2/m,2/m). 
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B. Kazakov Migdal Model 



The situation for the Kazakov-Migdal Model is significantly more complicated mostly 
due to the presence of the gauge field. One option for solving this model is to begin with 
the Functional Integral (0) and to define the two quantities E\ and Ga^ as follows: Let a 
and b be two adjacent sites on the lattice and let Uab be the gauge field on the link joining 
them then define the average < Q > of any quantity Q as 

<Q>= ^ / Y[d<p{x) n [dU{x,y)] Q exp [-^km] (16) 



We then define 



<x,y> 



which is expected to be independent of the chosen site a (in the limit of large volume) and 

which is expected to be independent of the chosen link (a, 6). It is useful to note that 
asymptotically, for large A, 

G,,^^ + ... (20) 



The next step is to write two equations analogous to Eq. one for E\ and one for Gx^ 



<x,y> '^^K'^/ij I -ry^j/ij 



and 



where S{(j), U) is the Kazakov-Migdal Action ([l|). 



Recall that in the limit of infinite N the functional integral is dominated by a single, 
translationally invariant Master Field Due to the Gauge Invariance of the Action we can 
choose $ to be diagonal with eigenvalues $j {i = 1..N) without any loss of generality. 

Before proceeding it is useful to define a quantity A{(f)) as follows: Consider first the 
quantity 

C^J{^) = I [rft/]|t/.,f exp {-Skm{U, $)) (23) 

Then define 

A,(<l>) = a,(<l>)<l>, (24) 



With this definition in hand one can derive the following equations (see Ref. pO[ for 
details) corresponding to Eqs. |2l| and ^ respectively: 



,2 ,l^„r($)-2DA($), 



Ei - (— Tr^ ^) = (25) 



(E. + - - (lT..IM?_imt,,,_l_^,.> . (26) 



Note A(0) is defined so that (■ ■ ■Uxy(p{y)U~y) is replaced by (■ ■ ■A($)) in the large N 
limit provided the full expression is gauge invariant and that Uxy does not appear in the 
expression again. At this stage we must make an important comment about notation. In 



Equation (^6]) the only average remaining is that over the Gauge fields which is done in the 
background field $. After doing this average all that remains is to perform the trace which 
can be thought of as an average over the distribution of eigenvalues of $. It is thus common 
to write, for example, V'{(f)) or A(</)) instead of V^'($) or Aj($) where we identify the label 
i with the eigenvalue so that $i = (p. In the case when the potential is the trace of a 
function of the matrix scalar field, V'{(j)) will have the same form as the derivative of the 
potential with respect to the matrix valued field. 

To proceed further we need to make a choice for the potential V{(f)). The solution to the 
above equations is simplest in the case of a Gaussian potential which is discussed in detail 
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in Ref. ||2^. In this case it turns out that A(0) is simply proportional to 0. As a result 
one obtains a quadratic equation for Ex whose solution (see |2y] and yields a semicircle 
distribution of eigenvalues where the edge of the distribution 2/M(m) is a calculable function 
of the parameter m in the potential. 

In this paper we shall discuss the more complicated case of a "Penner" potential. Previous 
work on this subject can be found in Ref. The case which can be solved involves choosing 
a potential for which 

V'{<P) - {2D - 1)A(0) = ^ + 5 (27) 

which has a pole at = ^ with residue q and an asymptotic value of i? at ^ oo. The idea 
is to use this form of V — {2D — 1)A(0) and to solve for the potential V which leads to this 
function. 

Using Eqs. ( ^5]) and (^) (see pOf for more details) we can derive the following equation 
which relates Gx^ to E\ and (which thus involves G at the singularity of V): 

_ (A - i)Ex - qG^^ 

- {X-0{Ex + ,-B)-q ^^'^ 
This can be written in a more useful form by using the asymptotic condition for Gx^- If we 
expand Gxu in Eq. (^8|) in large A using Eq. (^) we find that 

qG^^ = 1 + {B- fi)E^, (29) 

Substituting this back into Eq. (EBf) we relate Gxu to Ex and E^: 



_ {X-^)Ex + {fi-B)E,-l 



The next step is to notice from Eq. (|TBp that Gxu is symmetric under interchange of A 
and u so that Gxu = G,yx- Applying this condition to Eq. (|5D|) one finds, after some algebra, 
the following quadratic equation for Ex: 

(A - B)X -OEl + [(A -B){X- e)(e -B) + q{B - e) - (A - 0] Ex 

+ X{B-^) = constant (31) 



where the constant is independent of A. It is possible to determine this constant in terms of 
the mean value of 0; = / p{(p)(pd(f) using the asymptotic condition (|T9|). The equation for 
Ex then becomes: 

(A -B){X- El + [(A -B){X- 0(e -B) + q{B - - (A - 0] 



+ 



{X + ^){B-0 + ie-B') = (32) 



In principle it is necessary to determine the value of self consistently by extracting the 
density of eigenvalues p, which will depend on 0, from Eq. (|32D and then demanding that 
/ d(j)p{(f))(f) = (f). We shall see however that in many cases a properly normalized solution will 
exist for range of values of 0. The physical reason for this is that due to the singularities of 
the potential the eigenvalues (in large A^) can be arbitrarily distributed among two minima 
on different sides of a singularity. This point will be discussed further later in this paper. 

It is now possible to extract the functions A(0) and the potential V{(j)) by comparing the 
equation ( p2D with Equation (|25|). Notice that if V — 2DLambda would have only simple 
poles we would also obtain (from Eq. ( P5| ) a quadratic equation for Ex. By comparing these 
two equations one finds (see [^] for details) 

V\X)-2DA{X) = ^ + X^ + (^-0 (33) 



Using the ansatz (^71) we find 



A(A) = fr^ + e (34) 
which leads finally to the expression for V' 

^'<^> = + """^i^r" + (2« - + B (35) 

which results from our ansatz (pTj). This now yields the potential V for which we have found 
a solution. 

The preceding results thus establish an interesting relationship between the KM Model 
with two logarithmic singularities as in Eq. ( |55| ) and a ordinary one matrix model with the 
potential W{(f)) whose derivative is of the general form 
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(36) 



with 



ri + r2 



1; 



C 



(37) 



where ri is identified with q, rji with ^ and 772 with 5. A Matrix Model with the potential 
Vr(0) will have the same function E{\) and thus, in particular, the same distribution of 
eigenvalues as the KM Model. This One Matrix Model can be discussed for arbitrary values 
of ri,r2 and C (i.e. without requiring the conditions ( P?]) ) and it has many interesting 
aspects. For example it generalizes the Penner Model which contains a single logarithmic 
singularity to a model with two singularities and, as will be seen later, it admits a multi- 
phase solution space with non-trivial critical behavior. In Section |T| we thus study the 
double penner One Matrix Model whose results will later be applied to the KM Model. 



Before proceeding it is useful to recall that the infinite behavior of both the ordinary 
Matrix Model and the K-M Model can be described by the solution to an analogue mechan- 
ical problem. In the case of the ordinary Matrix Model the partition function Eq. (^ is 
written in terms of the eigenvalues 0j of the matrix 



In the large N limit in which the solution is given by the classical extremum of the action 
we see that we have an analogue mechanical problem of particles which are constrained 
to lie on a line at locations 0i ■ ■ ■ (pN- Each particle is subjected to an overall potential V{(f)i) 
and to a logarithmically repulsive two-body potential log — (pj)"^. 

The solution to the KM Model corresponds to another, more complicated, analogue 
mechanical model. If we look at the partition function in Eq. (^ and recall that in the 



C. Interpretation of Solutions 




N 



z 



i=l 




(38) 
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infinite N limit the (f){x) are independent of x and tliat on a square lattice each site has 
D independent nearest neighbors (where D is the dimensionality of the spacetime) then for 
infinite Eq-® is equivalent to solving the following one matrix model 



N 

1=1 



A2(0) exp 



-Nj2V{<f>. 



D 



N 

1=1 



exp 



-N E (Vi<p,) - D<p') -{D-1)Y: log (0. - <Pj)' + D log det,,e^(^-<^^)'/^ 



(39) 



This is once more a problem of particles on a line at locations (pi. Their central potential 
is now V{(j)) —Dcj)^ and their interaction is no longer a two-body interaction since it involves 
a determinant. We can, however, use the fact that the integral over the Gauge Group in 
Eq. (^ which lead to Eq. (^) is nonsingular when any two eigenvalues approach each other. 
Using this we see that the effective interaction is logarithmically repulsive at short distances 
(due to the extra factor of A^) and attractive at long distances if D > 1. 

The analogue mechanical problem presented in this section is a very useful tool for 
visualizing and checking the solution we obtain using the mathematical machinery of Matrix 
Models. 



III. THE DOUBLE PENNER MODEL 



A. Formal Solution 



In this section we investigate the large N solutions of the (non-KM) One Matrix Model 
with the potential (|36|) , for all possible values of the parameters ri, r2, t^i, r/2 and C. The 
basic equation is Eq. (|13D and leads to the quadratic equation 

El - W'{\)E^ + + = (40) 

X-r]i X-r]2 

with the Erf^ determined by the asymptotic expansion of Ex as 

riE^, + r2^r,2 = C 
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riEr,^r]i + r2-E^2?72 = C0 + ri + r2 - 1 (41) 



which allows the i?^. to be determined explicitly in terms of the mean value (j) of 
'nE^A I / (r?2 - 0)C + 1 - (ri + ra) 

The solution to Eq. ( ^OD is simply 



(42) 



Ex = l{ W'iX) 



\{W'{X)f-^j:^] (43) 
^ \ \ i=i ^ V' 

where 

W\X) = j:-^ + C (44) 

j=l ^ 'h 

is given by Eq. (|36D. Notice that we have chosen to write the solution to the quadratic 
with a "— " rather than with a "±" and to discuss the two solutions in terms of the possible 
branches of the square root. 

The correct choice of the branch of the square root in Eq. (H3) is very subtle, even in 



the case of the ordinary Penner Model (see [^). First of all we must choose the branch of 



the square root so that Ex satisfies the asymptotic condition (|T9|) : Ex ^ 1/ \ as \ oo. 
This will be satisfied provided the square root is chosen to have no branch cuts going out 
to infinity and that the square root approaches +C at infinity. There are however two 
other conditions which must also be satisfied. Notice from Eq. (^) than Ej^- is some finite 
number so that Ex is not singular at the singularities of the potential. We shall see that this 
condition is not automatically satisfied by the solution (^3]). Furthermore we would like the 
"Master Field" $ which is encoded by Ex to be a Hermetian Matrix with real eigenvalues. 
This requires the branch cuts to be on the real axis with a purely imaginary discontinuity 
so that the density of eigenvalues is real and positive. We shall see below that it is often 
impossible to satisfy all these conditions simultaneously though it may be possible to relax 
these conditions somewhat and still maintain an interesting solution. In fact we shall see 
that despite the fact that in the ordinary Double Penner model there are both one-cut and 
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two-cut solutions in the KM case in which the conditions (^) must be satisfied, only one- 
cut solutions will be possible (i.e. the second cut will correspond to two degenerate branch 
points leading to a zero density of eigenvalues). 

Before proceeding to analyze the branch cuts it is useful to write Eq. (^) in an alternate 
form. Note that the location of the branch points are found by finding the zeros of the 
function under the square root in Eq. (|43|) . This is a quartic equation for A. Let ^i, ^2, ^3, ^4 
be the solutions to this quartic equation. The requirement of a real, positive definite, density 
of eigenvalues will require these roots to be real though there may be some degeneracy among 
them. We thus expect, in general, both one-cut and two-cut solutions to the Double Penner 
Matrix Model. Let us assume, without loss of generality, that ^1 < ^2 < ^3 < ^4 Eq- (H ) 
for Ex can now be written 

As discussed above, the asymptotic condition on Ex requires the branch cuts of the square 
root to be chosen so that the square root approaches at infinity and that there are no 
branch cuts which go out to infinity. 

Implementation of the further conditions described above will be discussed in Sec. |III C 



Before doing so we should point out that the density of eigenvalues which is extracted from 
Eq. (|45| ) will always be normalized and have the correct value of provided only that Ex 
has no singularities and that there are no cuts going out to infinity. To see this recall that 
~ 1/A as A ^ 00. Thus 

^££.rfA = l (46) 

where C is a circle at infinity. If Ex has no additional singularities then this integral can 
be written as an integral over the discontinuities across the cuts of Ex which is just the 
total normalization / p{(f))d(f) which is thus equal to 1. The mean value of can similarly be 
calculated. From the asymptotic expansion (|19D of Ex we see that 

1 



^ . , XExdX = (f) (47) 
zm Jc 
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which is now vahd for the Ex in Eq. (^5]). Recall however that the parameters in 



depended on (f). The result (|47| ) is valid, however, independently of this value of (p. Again, 
if Ex has no additional singularities the contour can be deformed to a contour surrounding 
the cuts so that / (j)p{(j))d(j) = is this same value of which, as we recall, was extracted 
from the same asymptotic expansion of Ex- (If the cuts are chosen so that Ex is singular at 
either one or both of the rji then both the normalization and the mean value of will get 
an extra contribution from the singularity. It is not difficult to show, for example, that the 
normalization will get a contribution rj from a singularity at T]i. 

B. Review of Single Penner Case 

Before discussing the various possibilities in the Double Penner Case let us review briefly 
how the cuts work in the ordinary Single Penner Model. In this case the derivative of the 
potential is of the form 

W;i<P) = J^ + Q (48) 
and the quadratic equation for Ex resulting from Eq. (|13D has the solution: 



r , A' , Q 



+ Q (49) 



with 



rEr, = Q (50) 



The location ^± of the branch points occurs when the square root vanishes i.e. 

2-r 2 



When r > 1 there are no real solutions. This is the case when the potential Wp{(j)) has 
an attractive logarithmic singularity at t] with a strength greater than 1. When r < the 
potential has a repulsive singularity which results in a local minimum of Wp which occurs 

15 



at > ?7 if Q > and at < 77 if Q < 0. In this case it is easy to check that both ^± are on 
the same side of the singularity as this local minimum. In the case < r < 1 the potential 
has an attractive singularity and Wp has only a local maximum. In this case both branch 
points ^± lie side of the singularity opposite to this local maximum. 

In this Single Penner case we can, as discussed above, try to implement the three con- 
ditions: Ex 1/X a.s X 00; E\ nonsingular at rj; and p(0) real and positive with support 
on the real axis. The third condition implies that the branch points and branch cuts must 
lie on the real axis. This guarantees that p is real but not necessarily that it is positive. To 
analyze the other conditions we write, as in Eq. 



^^=2|a^ + ^-^^^ A^;^ j 

The asymptotic condition on Ex requires that there be no branch cuts at infinity and that 
the square root be positive for A > .^-t. This, together with the positivity requirement of 
p implies that the branch cut joins ^± in a straight line along the real axis. Furthermore 
the positivity of p requires that ^-t > 77 if Q > and < 77 if Q < 0. This is precisely 
the same result as Eq. (^Tj) so that in the single penner case we are guaranteed that the 
eigenvalue density will be positive if the branch cut joins the ^± in a straight line along the 



real axis. From Eq. (^91) we can determine the behavior of the square root in Eq. (|5^ at 
the singularity up to a possible sign. 

i^A-i^jr^li^r} asA-, (53) 

where the minus sign is used if the square root in ( ^2]) is positive which occurs ii t] > ^± 
and the plus sign is used if it is negative which occurs if rj < Thus we establish the 
following condition for E.^^ to be nonsingular: If 77 > ,^-1- then the sign of Q and r must be 
the same whereas ii rj < C,± the Q and r must have opposite sign. Thus, for the cancellation 
of singularities, we must have 

{v-^±)Qr > (54) 
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In the Single Penner Model the above condition is satisfied if and only if r < 0. This 
is evident from Eq. (|51| ) which implies that ^± > r/ if Q > whereas C,± < tj if Q < 0. 
Thus both cases require r < for the regularity of Ex at the singularity. This is physically 
reasonable since this is precisely the case in which the potential has a local minimum. 

The conventional way to analyze the case 1 > r > has been to relax the positivity 



condition on the density of eigenvalues (see pSf for example) and to allow the branch cut to 



go around the singularity so that the branch cuts goes from ^_ to ^+ by first going around 
7]. In the case 1 > r > this will lead to a distribution of eigenvalues with complex support 
but Ex will have both the correct asymptotics and it will be nonsingular at rj. One of the 
reasons this choice of cut is interesting is that it behaves nearly the same as a distribution 
of eigenvalues with a delta function singularity (of strength r) at rj and and the remainder 
of the eigenvalues along the real axis between the This can be seen by noting that if we 
compute the average of any analytic function /(A) 

I = [ dX p(A)/(A) (55) 

J cut 

then this integral will be independent of the precise path which the cut takes around the 
singularity. In fact by taking the cut from ^_ just below the real axis then around the 
singularity and back just above the real axis to ^+ one can easily check that 

/ = r/(r^)+ [^^ dX p{X)f{X) (56) 

with p normalized to 1 — r. This looks as if we could write 

p(A) = r6{X-v)+pW (57) 

with p having support on {^-,^+)- This is often called "condensation of eigenvalues" since 
some fraction of the eigenvalues "condense" at the singularity. Unfortunately this interpre- 
tation is not quite correct since if we were to evaluate the average of a nonanalytic function 
g{X) (for example if we were to evaluate the Free Energy which contains a logarithmic cut 
precisely in the region of interest) then the result would be dependent on the precise path of 
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the cut. In fact, for a general path, the Free Energy is not even real though some authors 
have chosen the contour in such a way that the Free Energy is real. 



C. Structure of Cuts in the Double Penner Case 

In the previous subsection we saw that even in the ordinary Penner model proper one-cut 
solutions which satisfy all our requirements do not necessarily exist and there are cases when 
either no solutions exist or when only very unusual solutions with "eigenvalue condensation" 



are present. We now continue with the discussion in Subsection (|III A| ) of the Double Penner 
Case where we shall find a similar situation. Let us focus attention on Eqs. (^31) and (^5|). 
As in the Single Penner case the requirement that Ex have no cut singularities at infinity and 
that p(A) be real and positive requires all the branch cuts to be on the real line. Recalling 
that C,i < C,2 < ^3 < ^4 we must then have one cut from to ^2 and the other from ,^3 to ^4. 
It is also clear from Eqs. p3| ) and p5| ) that if Ex is to be regular at the singularities they 
cannot be on a branch cut. Thus each singularity is either below ^1, between ^2 and .^3 or 
above ^4. Finally recall that the asymptotic condition on Ex requires that the square root 
in Eq. (|45|) must be positive for A > ^4. This in turn implies that it is negative between the 
branch cuts and positive for A < ^1. 

We are now ready to study the possible singularity of Ex near the r/j. Near the singularity 
at rji 

^^^Ij^ I T ^^^^r, I as A ^ (58) 

2X-rii [ \C\ n r]i-r]2 J 

(Near the singularity at ri2 simply replace 1^2 everywhere in Eq. (pSl).) The correct 
sign in the above equation depends on the sign of the square root near the singularity. The 
minus sign is to be used when the square root is positive i.e if rji > ^4 or r/i < ,^1 (in this 
case we shall say that the singularity is "outside" the cuts) and the plus sign must be used 
when it is negative i.e. if ^2 < ''7i < ^3 (in which case we shall say that the singularity is 
"inside" (i.e. between) the cuts). If the minus sign is used the singularity cancels provided 
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iVi ^ V2)Cri > whereas if the plus sign is used the cancellation occurs if and only if 
(^1 — V2)Cri < 0. Let us assume without loss of generality that 772 > Vi- (It is obvious that 
this can be done for the general double Penner case but for the KM Model the requirement 
that C > makes this not so obvious. There are however symmetries which relate the case 
T] > B to the case r] < B so that the assumption 772 > rji is completely general.) 
Our conclusion is then that: 

Cn > r]i inside < r]i < 6) -) 
Cri < rji outside (771 < ^1 or ^4 < r^i) +) 
Cr2 > 1J2 outside (772 < ^ or ^4 < 772) (^/ +) 
Cr2 < r]2 inside (6 < m < 6) -) 
(59) 

(If ?7i > 7/2 this would of course be reversed.) We shall consider the case C > throughout. 

Although the above conditions guarantee that Ex is nonsingular it does not guarantee 
that p is positive. The best way to see this is to look at Eq. (^Sj) . It is clear that if there 
are two cuts with no singularity separating them or with both singularities separating them, 
then one will have a positive p and the other a negative p. If we define the regions I, II and 
III as the regions A < 771, 7/1 < A < 7/2 and A > 772 respectively then one cut must be in 
region II and the other must be in region III. There is however a loophole. If one of the cuts 
is degenerate (i.e. ^1 = ^2 or ^3 = ^4) then p is zero along the degenerate cut and it does not 
matter if it is in the wrong region. We can thus have single cut solutions as follows: Either 
the cut (^1, ^2) is in region II or the cut (,^3, 1^4) is in either region I or III with the other cut 
being degenerate. 

We now proceed to discuss the various cases individually. 

D. The Case ri < and r2 < 

The simplest situation occurs when both ri < and r2 < 0. This corresponds to a 
potential W{(j)) which has two repulsive logarithmic singularities and two local minima. One 
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minimum is between the two singularities and the other is above r]2 (since we have assumed 
that C > 0).(We call this the "two up" or "2u" potential.) In this case the condition ( ^9]) for 
a nonsingular Ex is that the branch cuts lie in the same region as the minima of the potential 
namely in regions II and II. This is what we expect since the eigenvalues of a Matrix Model 
are expected to be distributed about the minima of the potential as discussed in Sec. [11 C 



The only remaining question is whether the branch cuts which are a solution to the quartic 



equation resulting from Eq. (^3|) do in fact lie in the correct place. It is possible to show 
that they always do which also guarantees us that the density of eigenvalues is everywhere 
positive. (Recall from the previous section that when positivity of a double cut solution 
requires the cuts to be in regions II and III when C > 0.) 

The solution to the quartic equation leading to the branch points ^1,^2,^3,^4 of Ex in 
Eq. ( ^3]) is in general quite complicated. The simplest way to see the types of solutions 
which are possible is by a graphical method. Our goal is to solve the equation 

^ ■ +cV-4f-^ + ^;^Uo (60) 



\A-?7i X-r]2 J \>^-Vi A-r^a 
with TiEi and ra-Ea given in Eq. ( ^Hj ) in terms of 0. Without loss of generality we may 
choose r/i = and call ri2 = rj. Let us also call riEi = 6 so that ra-Ea = C — 6. Eq. ( pDD is 
thus equivalent to: 

[ri (A - + r2 (A - r],) + C (A - r,,) (A - r/2)]' = 4((CA - 5ri) ((A - r],) (A - r/2) (61) 

The left hand side (LHS) has zeros precisely at the extrema of the potential W{X) which we 
call Ai and A2. In our present case (ri < 0, r2 < 0) we have < Ai < 77 < A2. We now sketch 
both the LHS and the RHS of Eq. (|6TD for various values of 6 (which, we recall, is related 
linearly to 0). The LHS is a quartic with the two degenerate roots Ai and A2 whereas the 
RHS is a cubic which goes like 4CA^ as A — > +00 and with roots at rji, 772 and Srj/C. We 
now notice the following: 

For 6 < there is always a pair of roots in region HI (A > rj) but there are never any 
roots in region II. Depending on the values of the various parameters there may be a pair 
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of roots in region I. These would however lead to a negative density of eigenvalues and, 
correspondingly, the p(A) for A > would have J^^^Q p{X)d\ > 1. In case there are no roots 
in region I the two additional roots are complex and, consequently, the density of eigenvalues 
for A > will not be normalized to 1. 

When 6 is slightly positive, the roots in region III persist but there are no roots in region 

I nor in region II. At some critical 6ci which lies in the interval (0, Ai) a pair of roots begins 
to appear in region II. This pair of roots persists for all 6 > 6)cl. At some critical 6c2 > A2 
the pair of roots in region III disappears. 

Thus for all values of (p for which 6ci < 6 < 6c2, Ex has two branch cuts in regions 

II and III for which the density of eigenvalues is positive, normalized to 1 and for which 
/ (j)p{(f))d(f) = as required. For the two special cases 6 = 6ci and 6 = 6c2 two of the branch 
points become degenerate and the above two-cut solution reduces to a one-cut solution. 
Notice also that for all the above cases p(A) ~ |A — near any branch point ^. 

The "physical" reason for the existence of this large class of classical solutions can be 
seen by referring to the mechanical analogue problem discussed in Sec. |11 C\ Since there is 



an infinite barrier separating the regions II and III we expect that a solution will exist with 
any number rii of particles in region II and N — rii particles in region III. The degenerate 
cases 5 = 6ci and 6 = 6c2 correspond to the cases when all particles are either in region I 
or in region II. 

E. The Case ri > and r2 > 

The next case we consider is the case in which both ri and r2 are positive. (This case 



could correspond to a KM Model provided the conditions (|37D were satisfied.) 

The next case we consider is the case in which both ri and r2 are positive. (This case 
could correspond to a KM Model provided the conditions in Eq. ( P7| ) were satisfied.) In this 
case the potential has no minimum but it has two maxima at points Ai and A2 in regions I 
and II respectively. We expect no real normalizable solutions in this case though there may 
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be cases in which "eigenvalue condensation" occurs as was discussed in Sec. [Ill B| for the 
Single Penner case. 

Using the conditions given in Eq. ( pQ]) (and recalling that C > 0) we see that in order 
to avoid singularities of Ex at rji and 772 we must have one cut in region I and the other in 
region II (i.e. in the same regions as the extrema of W{X)). This, unfortunately, leads to a 
negative density of eigenvalues in both regions I and II. We thus conclude that there are no 
normalizable one-cut or two-cut solutions in this case. 

In order to examine the possibility of "eigenvalue condensation" we use the graphical 
approach discussed in Sec. [Ill D| . As in this previous section we assume without loss of 
generality that 771 = and we call ri2 = rj and we conclude as follows: For many values of the 
parameters there will not be four real branch points. In case four real branch points exist 
there are two possibilities. Either the two branch cuts are in regions II and III or in regions 
I and III. Let us begin by discussing the case when the cuts are in regions II and III. In this 
case p(A) is positive along both cuts but Ex is singular at both rji and r/2 unless we deform 
both branch cuts in regions II and III to surround the singularities rji and 772 respectively. 
Recall from the single Penner case that for the purposes of computing averages of analytic 
functions we can choose a contour which circles around each singularity but otherwise goes 
along the real axis between (^1 and ^2) and between (^3 and ^4). The singularities at 771 and 
ri2 contribute an amount ri and r2 respectively to the normalization. We thus expect that 
such solutions should exist only if ri + r2 < 1. 

Finally we examine the possible case in which the cuts lie in regions I and III. In this 
case p(A) is positive in region III but negative in region I. Thus, in order to have a sensible 
solution, the two branch points in region I must be degenerate (^1 = ^2) so that the negative 
p in region I is of no concern. We will however have a singularity of Ex at 772 unless the cut 
in region III surrounds the singularity 772- In this case the singularity will contribute r2 to 
the normalization so we expect such solutions only if r2 < 1. 
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F. The Case ri > and r2 < 



In the case ri > r2 < the potential has a minimum in region III and a maximum in 



in regions I and III. In this case the density of eigenvalues p(A) will be positive in region 
III but negative in region I. It thus follows that the branch points in region I must coincide 
{^1 = ^2)- It is easy to see using the graphical method described in the previous sections 
that, independent of the parameters of the potential (provided ri > r2 < ) such a 
solution to our quartic equation always exists. This leads in every instance to a normalizable 
distribution of eigenvalues in region III for which the eigenvalues are distributed about the 
minimum of the potential. 

We now consider the more general possibility in which we allow "eigenvalue condensation" 
by relaxing the condition (|59|). Depending on the values of the parameters in the potential, 
there may be no real set of branch points or else the branch cuts will be in regions II and 
III. This leads to a positive density of eigenvalues along both cuts but Ex will be singular 
at rji unless the branch cut in region II surrounds 771. Since this singularity will contribute 
ri to the normalization, we only expect this solution to occur if ri < 1. 



region I. The condition for the existence of a nonsingular solution is that the cuts lie 



G. The Case ri < and r2 > 



Finally we turn to the case ri < r2 > 0. In this case there are four possibilities: 




4|ri|Cr7 > (ri + r2 — Cr/)^ 



(62) 



In this case W{X) has no extrema at all. 

In the remaining cases 4|ri|C?7 < (ri + r2 — Crff' . 



Ti + r2 > Crj 



(63) 



W{X) has first a maximum then a minimum in region I. 



(c) 



— Ct] < ri + r2 < Crj 



(64) 
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W{X) has first a minimum then a maximum in region II. 

(d) ri + r2 < -Cr] (65) 

W{X) has first a maximum then a minimum in region III. 

In the case (a) we expect no nonsingular solution. Whereas in the cases (b),(c) and 
(d) we guess that if the minimum is not too shallow we should get a normalizable solution 
distributed about this minimum. The condition ([59|) for nonsingular solutions in this case 
allows the branch cuts to be either both in region I or both in region II or both in region 
III. 

We begin by analyzing case (a) in which Vr(A) has no extrema. In this case, the graphical 
method implies that the four real solutions to the quartic equations, if they exist, may be 
either in regions I and III or in regions II and III. It can further be shown that they cannot 
be in regions I and III. Thus they can only be in regions II and III. We thus conclude, as 
expected, that in this case we do not have a nonsingular solution. However in case there are 
two pairs of roots in regions II and III, the resulting branch cuts lead to a real eigenvalue 
distribution for which Ex will be singular and 772 unless the cut in region III is deformed to 
circle around ri2. This will lead to "eigenvalue condensation" at ri2 which can only occur if 
r2 < 1. 

Next we look at case (b) for which W{X) has first a maximum then a minimum in region I. 
In this case the quartic equation has four roots in region I (for suitable values of 6) provided 
only that (ri + r2 — C?7)^ + 4riC?7 is not too small (i.e. the minimum of is not too shallow) . 
The cut connecting ^1 and ^2 must be degenerate since it would lead to a negative p{X). 
We thus have a normalizable, nonsingular one-cut solution which can be seen (using, for 
example, the graphical method) to lie around the minimum of W. A precisely analogous 
situation occurs for cases (c) and (d). We have a normalizable, nonsingular distribution of 
eigenvalues near the minimum of the potential provided the minimum is not too shallow. 
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IV. THE KAZAKOV MIGDAL MODEL 



We saw in the previous section that the ordinary "Two-Pole" Penner Model has a rich 
variety of one-cut, two-cut and singular solutions. In this section we apply these results to 
the KM Model. According to the results of Sec. [11 B| the solution to the KM Penner Model 



with the potential ^(A) given in Eq. (^) is related to an ordinary (non KMM) "Two-Pole" 
Penner Model with W{X) given in Eq. ( PU]) provided (P^) 

ri + r2 = l; V2 - Vi = C (66) 

(Recall that ri = q, rji = ^ and 772 = B.) For simplicity we shall call ri = r which is, of 
course, also equal to q. The results of Sec. |T| can now be applied directly to the KM Model. 

The most obvious conclusion we can draw immediately is that the case ri < r2 < 
is not applicable to the KM Model. This was the only case in the previous section which 
admitted a nonsingular two-cut solution. We thus conclude that in the KM Double-Penner 
Model we have at best only one-cut solutions. 

Consider first the case r > 1. In this case r = ri > 1 > and r2 = 1 — r < 0. This case 
was discussed in Sec. [Ill h] . For all such values of r and for all C = rj = ri2 — rji we have a 
normalized one-cut solution in region III. In fact the degenerate branch points C,i = ^2 < 
and ri2 < ^3 < X2 < ^4 where A2 is the location of the minimum of W{X). The alternate 
solution in which there is condensation of eigenvalues at the pole rji cannot occur in the KM 
case since r = ri > 1. 

It is instructive to examine the KM Potential in this case and then to compare our results 



with the physical expectation described in Sec. |II C[ . From Sec. [II C[ we recall that we are 
interested not simply in the potential V which appears in the KM action but in V{X) — DX^ 
(see Eq. (^)). Choosing, without loss of generality, ^ = and calling B = r] we have 

V^'(A) -2DX = ^+ (2^- l)(^-l) ^^_2DX (67) 
A A — Tj 



V{X) - DX^ = rlog(A) + {2D - l)(r - l)log(A - r]) + r]X - DX^ (68) 
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Noting that 



7]\ - = -d(\- + constant (69) 



2 



we see that, unhke W{X) which has a maximum in region III, V{X) — DX^ has a maximum in 
this region. This is an example of the situation discussed in Sec. [11 C| in which the eigenvalue 
distribution has its support at the maximum of a potential. The eigenvalues can "straddle" 
the maximum due to the long range attraction of the eigenvalues which is a result of the 
integral over the Gauge Fields. 

The next case we consider is when r < 0. In this case ri = r < and r2 = 1 — r > 1. 



This situation was discussed in Sec. Ill G. First note that there can be no eigenvalue 



condensation in this case since r2 > 1. Furthermore Eq. (|62|) implies that if 

there are no real normalizable solutions. In fact even if |r| is somewhat less than this limit, 
the minimum of the potential W{X) is too shallow to admit normalized solutions. When \r\ 
is sufficiently small we can look for solutions by first noting the location of the extrema of 
W{X) using Eqs. ( p5|J5^ , p^D . If 77 < 1 the requirements of case (b) of Sec. [ill G| is satisfied 



so W{X) has a minimum in region I. It can be shown that in this case there are never four 
real roots in region I. (One way to see that this has to be the case is by noting that if such a 
solution did exist then (p would have to be in region I. This would imply that 6 = rj — (f) > 1]. 
The graphical method of the previous section then shows that there can be no roots in region 

I-) 

If > 1 we satisfy the requirements of case (c) resulting in a density of eigenvalues whose 
support is in region II. There is however a critical value rj^r > 1 of rj with the property that 
for rj < rjcT there is never a normalized distribution of eigenvalues (i.e. one never has two 
cuts in region II one of which is degenerate). For rj > rjcr a proper solution exists provided 
|r| < |rc(?7)| < {1] — l/riY/A. The curve rc{ri) can be computed with the aid of Eq. (^Sf). 
Using the fact that E\ is nonsingular at both A = and ai X = r] and using the fact that 
Ex 1 / X a,s X ^ 00 one derives the following equation: 
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(71) 

The critical line will occur when three of the are degenerate. This leaves three equations 
with only two .^'s unknown and allows us to determine r in terms of 1] with the result that 

When r < rc{ri) the KM potential V^(A) — DX^ has a minimum in the support of p. Notice 
that there is no case when r < for which the minimum of W is in region III. 

Finally we consider the case < r < 1. In this case 1 > r = ri > and 1>1— r = r2>0 



which is discussed in Sec. [Ill E| . There are no nonsingular solutions in this case. Unlike the 
previous examples this is a case in which the ordinary Penner model can have "eigenvalue 
condensation" since both ri and r2 are < 1. In fact we found cases in which "condensation" 
occurred just aX X = rj and cases in which is occured at both A = and A = 77. Unfortunately 
the KM Potential V - DX^ has a maximum at A = 77 and thus eigenvalue condensation 
cannot occur at this point. The only possible explanation is that the procedure which 
lead from the Ordinary Penner Model to the KM Penner model does not work when the 
eigenvalue distribution is singular. This is not surprising, especially in light of the problems 
with interpreting these as real distributions. We thus conclude that for the case < r < 1 
there are no solutions at all. 

In summary we see that we never have eigenvalue condensation in the KM Penner Model 
and that nonsingular normalized solutions exist for all rj when r > 1 and for |r| < |rc(r7)| < 
{rj — l/riY/A when r < and rj > 1. The phase diagram for this model is shown in Figure 1. 
There are several interesting critical lines. Along the line 1] = the potential W'{X) = 1/A. 
This is an single-pole Penner potential with a critical value of the coupling (r = 1) but 
without a linear term. The line r = 1 [i] > 0) yields a potential W'{X) = 1/A + 77. This 
is again a critical single-pole Penner potential but this time with a linear term. The line 
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r = is also a critical Penner potential though centered at X = rj and with a linear term. 
Finally there is the critical line r = rc{ri) given by Eq. ( [71]) at which the potential in region 
II become sufficiently deep to admit normalized solutions. 

The behavior of physical quantities such as the susceptibility 

X = -^^nr) (73) 

where F is the Free Energy near the various critical points or lines is a subject of much 
interest in Matrix Models. In the case of the KM Penner Model this was studied in great 



detail by Makeenko ||2T| who computed the susceptibility as well as the various critical 



exponents of the model. 
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V. DISCUSSION 



In this paper we have studied in detail the large solutions to a Kazakov-Migdal Model 
with two logarithmic singularities 

V{(f)) = r log0 + {2D - l)(r - 1) log (p ~ r] + r] (74) 

and the solutions to the related problem of an ordinary Matrix Model with the Double 
Penner Potential 

W{(j)) = nlogcj) + r2\og(j) - 7] + C (75) 

The ordinary Double Penner Model has a rich phase structure which includes regions in 
parameter space in which there are one-cut solutions and two-cut solutions, regions in 
which there are no solutions and regions in which the solutions are singular with "eigenvalue 
condensation" at the poles. The KM Penner Model on the other hand has either one-cut 
solutions or no solutions at all. Its phase diagram is shown in Fig. 1. 

The method employed in this paper, which was first used for the KM Model in Ref. |P0|| , 



emphasizes the power of the Path Integral in solving this difficult mathematical problem 
which otherwise would have involved finding the extrema of Actions involving the logarithm 
of the Itzykson-Zuber determinant in Eq. (^). The techniques of Ref. reduce the 
problem to solving an ordinary quartic equation. 
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FIGURE CAPTIONS 
Figure 1: Phase Diagram for the Kazakov-Migdal Penner Model 



31 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9501037vl 





SOLUTION IN REGION III 


r=l 










NO SOLUTIONS 




r=0 




NO MINIMA 






NO 


SOLUTib?^ \v 


SOLUTION 




MINIMUM 
IN 

REGION I 


^\ MINIMUM \. 

\ IN \ 
NO MINIMUM \ ^(.joj^ 

\ 


IN 

\ REGION II 



r|=0 r|=l 



Figure 1 



